—_ 


CONTENTS 


Abstract 

Problem Status 

Authorization 

INTRODUCTION 

COMPUTATION PROCE DURE 

METHOD OF CHECKING 

RANGE OF TABULATION 

COMPARISON WITH OTHER PUBLICATIONS 
REFERENCES 

TABLES OF EIGENVALUES OF THE EXTERIOR 


ACOUSTIC RADIATION PROBLEM FOR 
PROLATE SPHEROIDAL COORDINATES 


ABSTRACT 


The wave equation in prolate spheroidal coordinates was 
separated into radial and angle functions. The differential equa- 
tion satisfied by the angle functions was written in the form of an 
eigenvalue problem, that is, as a linear operator operating on 
eigenfunctions to yield the same eigenfunctions multiplied by cor- 
responding eigenvalues. The eigenvalues were numerically cal- 
culated by use of Galerkin's method. This method reduces to the 
evaluation of the characteristic roots of a large matrix. An 
80 by 80 matrix is chosen and a detailed calculation on a high- 
speed computer leads to a tabulation. The table of prolate eigen- 
values published here has the range m = 0, 1, 2; = m(1) m+ 49; 
h = 0.1 (0,1) 0.9, 1.0 (0.2) 8.0, 10.0, 20.0 (20.0) 100.0. The pre- 
cision is 21 significant figures. 


PROBLEM STATUS 
This is an interim report on one phase of the problem; work 
on other phases continues. 
AUTHORIZATION 


NRL Problem S01-28 
Project RR 102-08-41-5225 
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TABLES OF EIGENVALUES OF THE WAVE EQUATION 
IN PROLATE SPHEROIDAL COORDINATES 


INTRODUCTION 


In the solution of Helmholtz's wave equation (\?+ k)w- 0 by separation of variables 
in spheroidal coordinates (¢.1.0) the following equations occur: 
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Here A is a Separation constant, h kd 2, k .,, ¢ is the interfocal distance of the 


spheroid, F is a radial function, 5 is an angle function, and # is a rotational function. 
Equations (1) and (2) have the form of the Liouville equation 
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Since Eqs. (1) and (2) are identical, except for range of the variables , and ., we con- 
sider only Eq. (2). Here we seek the real eigenvalues A 4, () = m.m+1....) for which 
the solutions of Eq. (2) over a range are nontrivial and nonsingular. The singularities of 
Eq. (4) are known to occur at p(7) 0, ie., at » ‘1, The range of interest is thus 
-1< = +1, or parametrically ,» cos 4, 02 -:. The boundary conditions are of 
the conventional type, ata,’ Oat y=1, and pyr fy’ Oat , -1. For these 


boundary conditions there will be a lowest eigenvalue 4, which will be smaller in mag- 
nitude than any other 4,,, and all other eigenvalues can be ordered to form a monoton- 

ically increasing sequence 4,,;,,, ° 4,,- To find the eigenvalue 4, we use the method 
of Galerkin (1); i.e., we assume that the corresponding eigenfunction can be expanded in 

a finite orthogonal series {.,(,)} over the range -1-  ~ 1 and write 
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in which , is a solution of Eq. (4) and associated homogeneous boundary conditions in 
the range -1_< , ~ +1 with known eigenvalue \,. Substituting Eq. (5) into Eq. (2), multi- 
plying by , and integrating over the range according to Galerkin's procedure, one gets 
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An appropriate choice of ¢, for the range of , 18 the associated Legendre function 
P™ (1) Which satisfies Eq. (2) for the condition h=- 0. This condition requires that 
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Substituting 4, = P,”,, and 4, - P,” in Eq. (6) and using Eq. (7), it is easily seen that 
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Eq. (6) can be written in the form 
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where 
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The quantity L can be regarded as a large (though finite) matrix having elements L,,,, 
with eigenvalues ,. In this light, Eq. (8a) is an eigenvalue problem and the entities C,, 
are the components of an eigenvector corresponding to each eigenvalue a, Now let 
m+rz=p for convenience and write the normalization and recursion relations of P,” in the 


forms 


ste ate 1. p=s, 
| pee BS” inde (9a) 

BS : 0. pts, 

% P m ' 
pm p - hee, 2 2 ptm)! — Gb) 
Paes YO = Bert Corey! ” 
(2p 61) 9 Pm Feo oT es Cacia Rs (9c) 
In normalized form the last equation reads 
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The products 4(q)B(p) may have the following two forms depending on the procedure of 
normalization: 
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or 
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Equation (9f) is designated as the nonnormalized form, and Eq. (9g) is the normalized 
form. 


Expanding Eq. (8b) by use of Eq. (9d) one obtains further 
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It is Seen that any three consecutive terms of Eqs. (8a) satisfy the recursion relation 
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This is the nonnormalized recursion relation between the constants C_,. 
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We return to Eq. (9d), which is in the form 
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Substituting this formula into Eq. (8b) one obtains 
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Using the properties of normalization shown in Eq. (9g) this reduces to the recursion 
relation 
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This is the normalized form of the recursion relation between the constants C,,. 


In the computation of eigenvalues it is advantageous to use a matrix formulation 
with normalized elements as given by Eq. (12). However, in the determination of eigen- 
vectors whose components are the constants C,, it is convenient to obtain these in un- 
normalized form by use of Eq. (11) with appropriate choices of ,, and then to normalize 
them to suit a particular requirement. In this report the following normalization for the 
components of the eigenvectors is used: 
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This choice of normal zation (see Ref. 2) is arbitrary but very suitable to numerical 
computation. 


COMPUTATION PROCEDURE 


All elements in matrix L of Eq. (8b) are real. From Eqs. (11) and (12) it is seen 
that the subscripted elements s=p, s = p+ 2, and s = p- 2 are non zero ile all 
other elements of the matrix vanish. Furthermore, one notes that the matrix is Sym-~- 
metrical; i.e., C, , = C, , Based on these observations, therefore, we are required to 
find the eigenvalues , of a real symmetric matrix. To assist in the numerical compu- 
tation we have used the standard IBM Share-Distributed 704-709 Fortran program 664 
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by Burton S. Garbow of Argonne National Laboratory (CO-OP ID F4 UCSD 1 EIGEN). 
This subroutine finds all scalar solutions to the matrix equation 4x = ax, where Aisa 
real symmetric matrix. Although the eigenvalues of the separated Helmholtz equation 
are infinite in number, finite sized matrices yield excellent and well converged eigen- 
values provided the selected partition of the infinite matrix is "large enough". In our 
application we use an 80 by 80 matrix to compute 50 eigenvalues in double precision (21 


significant digits). 


METHOD OF CHECKING 


In Eq. (11) we designate the ratio C,,C, ,.. by the symbol ",". Inthe limit r + ~ 
we require that '/" vanish. Changing r to r+ 2 so that the ratio of constants is #7 , 
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one sees by the recursion relation Eq. (11) that #7 , can be written as an infinite con- 
tinued fraction [ i by the permutation of subscripts, i.e., increasing subscripts so 


r+2 
that the ratios read C, .,4C, p42) Cs.s., ©s c.g etc. Another set of permutations is 
the set of decreasing subscripts C,, C. __,, C, ,., C, -.4, etc. This latter set is lim- 
ited by the requirement thatC.,- 0 if : » 0. We canthus write "as a finite con- 


tinued fraction ["” ,] , using this latter set, i.e., the decreasing subscript set. Clearly 
one requires that 
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This is a transcendental equation in the unknown , whose roots are the eigenvalues 4... 
A skillful procedure for finding these roots has been devised by Bouwkamp (3). He takes 

the first variation of U(,) due to a variation 5,, devises an equation for 5, in terms of 

the 7," ratios, and by iteration of Eq. (14) using \15\,:4,,, etc., finally causes to 

converge to the true value 4.. This method however requires a good first guess on a ,, 

which can then be used as a basis for iteration. In the event that this guess is poor, the 

method yields a different eigenvalue \,,,, 7 = 2,4,6,... . Comparison of the matrix 

method noted above with the method of Bouwkamp shows that the matrix method agrees 

with Bouwkamp's method to all 21 significant digits for h << 10. For h > 20 Bouwkamp's 

method yields eigenvalues but not in a monotone series. Because these eigenvalues are 

sometimes identical ior different ™, !, h, it is seen that Bouwkamp's method fails for 

large h. The failure is attributable to wrongly chosen starting values. To overcome 

this difficulty an 80 by 80 matrix is used to obtain eigenvalues which provide starting a 
values for Bouwkamp's method. This procedure of combining the two method: yields 50 

correct eigenvalues up to 4 = 80. For h = 100 the same method yields approximately 30 

eigenvalues with 21 significant digits. This limit in available eigenvalues with required 

accuracy is due to limitation of the Share program, 


RANGE OF TABULATION 


The range covered by these eigenvalue tables is as follows: 


m= 0,1, 2 
i = m (1) m+ 49 
h = 0.1 (0.1) 0.9, 1.0 (0.2) 8.0, 10.0, 20.0, 40.0, 60.0, 80.0, 100.0 


The accuracy is considered to be 21 significant figures. 


6 HANISH AND KING 


COMPARISON WITH OTHER PUBLICATIONS 


The existing literature contains the following ranges of tabulation of prolate sphe- 
roidal eigenvalues. 


Chang and Yeh (4): 
m = 0 (1) 9, ¢ =m (1) 9, 4 = 0.1 (0.1) 10 (1.0) 100 
Accuracy = 12 significant figures. 
Stukey and Layton (5): 
m= 0 (1) 9, ¢ = m (1) m+9, h = 0.25 (0.25) 10.0 (1.0) 20.0 
Accuracy = 11 significant figures, 


The number of values of ¢ (i.e., the number of eigenfunctions for fixed m) in these 
latter collections is too limited to allow complete field calculations based on orthogonal 
expansions in eigenfunctions. For example, when m- 0, the Chang and Yeh tables contain 
ten values of ¢. Since in many physical problems only even or odd integers are per- 
mitted, this number has been found to be too small to permit convergence, except for 
very small values of h. The NRL tables in contrast provide 50 values of ¢ for has 
large as 80, 


The increase in accuracy (from 12 to 21 figures) in these NRL tables is motivated 
by a desire to permit accurate determination of the components of the eigenvectors (that 
is the expansion constants d,(h m?)), by a three-term recursion formula (6), Low 
accuracy of the eigenvalue results in poor d, constants. By increasing the accuracy of the 
eigenvalue to 21 places, we are able to obtain well over 100 «,, constants having accuracy 
comparable to that of the eigenvalue. Accurate determination of functions involved in 
field calculations is not possible without this capability. 
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TABLES OF EIGENVALUES OF THE 
EXTERIOR ACOUSTIC RADIATION PROBLEM 
FOR PROLATE SPHEROIDAL COORDINATES 


Range of Parameters 


0.1 (0.1) 1.0, 1.2 (0.2) 8.0, 10, 20 (20.0) 10.0 
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